Generalized scaling relations for unidirectionally coupled nonequilibrium systems 
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Unidirectionally coupled systems which exhibit phase transitions into an absorbing state are 
investigated at the multicritical point. We find that for initial conditions with isolated particles, each 
hierarchy level exhibits an inhomogeneous active region, coupled and uncoupled respectively. The 
particle number of each level increases algebraically in time as N{t) ~ with different exponents rj 
in each domain. This inhomogeneity is a quite general feature of unidirectionally coupled systems 
and leads to two hyperscaling relations between dynamic and static critical exponents. Using the 
contact process and the branching-annihilating random walk with two offsprings, which belong to 
the DP and PC classes respectively, we numerically confirm the scaling relations. 
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Nonequilibrium absorbing phase transitions (APT) 
have been studied extensively over the last decade in the 
hope of understanding better nonequilibrium continuous 
phase transitions P, Q ■ Also APT are one of the sim- 
plest and natural extensions of the well-established equi- 
librium transitions to nonequilibrium systems. In APT, 
a system undergoes a continuous phase transition from 
an active into an absorbing state where no dynamics oc- 
cur furthermore. Models exhibiting APT describe wide 
range of phenomena such as epidemic spreading, catalytic 
chemical reactions, surface growth, wetting and roughen- 
ing transitions, self-organized criticality, and transport in 
disordered media 0, U • 

The basis of theoretical and numerical analysis is the 
understanding that the concept of universality plays a 
central role in classifying nonequilibrium phase transi- 
tions as it does in thermal equilibrium. However, surpris- 
ingly few universality classes for APT have been identi- 
fied so far, and - as a unifying theoretical framework is 
not available - we are still far from a systematic classifi- 
cation of APT even in one spatial dimension. The best- 
studied directed percolation (DP) class includes the sys- 
tems which have no symmetry between absorbing states 
and no conservation laws in the order parameter 0, |3| ■ 
Another well-known class is the parity-conserving (PC) 
class which describes models with Z2 symmetry between 
absorbing states or parity conservation of the number of 
order parameter |3,y|- Recent studies show that the con- 
servation of total number of particles and couplingwith 
auxiliary diffusion fields also form distinct classes!!, , 
or give rise to unexpected new critical exponents 

Several critical exponents characterize the off-critical 
and critical behavior of APT. The exponents (3, 
characterize the scaling behavior of steady-state particle 
density p'', the correlation length ^, characteristic time r 
near criticality in thermodynamic limit as ||3| 



(1) 



probability of a single particle Pg (t) , density of particles 
p{t), the number of particle N(t) averaged over all sam- 
ples and spreading distance R{t) averaged over survival 
samples as 15] 



Ps{t) ^ t-\ p{t) ~ i-", N{t) ^ t\ R{t) ^ t 



l/z 



(2) 



with a — I , z = v^^/v^ and d is the spatial dimension. 
The well-known scaling relation between these exponents 

is Hi 



rj — d/z — a — S 



(3) 



Except in systems with infinitely many absorbing states, 
we have a = S 

Rather than exploring new systems exhibiting un- 
known critical behavior, coupled systems of known uni- 
versality classes have been studied recently as an another 
direction of searching for new universality classes. A 
coupled systems is a multi-species system in which each 
species is coupled to the others in certain ways. Among 
possible ways of coupling such as bidirectional, cyclic 
and unidirectional coupling in linear or quadratic way 
Hi [13 1 a recent study on a hierarchy of linearly unidi- 
rectionally coupled DP systems shows that a new critical 
behavior emerges at the multicritical point where criti- 
cality of all levels coincide . The multicritical behavior 
is characterized by a varying order parameter exponent 
/3, taking a different value in each level of the hierarchy, 
while the exponents i'± and ^| | are those of the DP class 
in all levels. Also rj and i5 vary according to the levels 
of the hierarchy, but these exponents do not satisfy the 
ordinary scaling relation of Eq. (3) in each higher level 
k{> 1). Instead they are suggested to satisfy ^ 



rik ^ d/zk- ttk - 5i 



(4) 



where A is the distance from criticality. At criticality, 5, 
a, 77 and z characterize the dynamic scaling of survival 



with ak = 5k- In Eq- (4), the exponent of Ps{t) of the 
first level, 5i replaces 5k of A;' level. The appearance of i5i 
in Eq. (4) is unusual and there has been no clear under- 
standing for that. It should be noted that Eq. (4) comes 
from taking an average where one uses configurations in 
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which each level survives regardless of the survival of the 
first level. We shall refer to this usual average as slave 
average to distinguish it from our average method dis- 
cussed below. 

In this paper, we report generalized scaling relations 
for linearly and unidirectionally coupled systems belong- 
ing to same or different universality classes. We propose 
an explanation for the appearance of in Eq. (4). We 
also show that Eq. (4) is not satisfied in general cou- 
pling systems because it is a consequence of the usual 
slave average which does not take into account an impor- 
tant inhomogeneity in the distribution of particles in the 
higher levels of the hierarchy. In addition, our scaling 
relations provides insight in the origin of markedly large 
corrections to scaling which make estimations of critical 
exponents difficult. 

We consider a hierarchy of an arbitrary number of cou- 
pled systems with coupling A B — > C - -. For the 
derivation of scaling relations, we are interested in the 
dynamic behavior of the fc'th level starting with a single 
particle on the first level. Then on each level, a clus- 
ter is created and spreads. Due to the coupling, the 
size of the cluster on each level is always larger than 
those of the lower levels. It means that the cluster of the 
fc'th level (with size Rk{t)) is divided into two parts, the 
coupled and the uncoupled region where the dynamics 
on level k evolve autonomously (Fig. 1). In the cou- 
pled region (with size Rc{t)), the (k — l)'th level feeds 
particles to the k'th level so that Rcif) is actually the 
spreading distance of the (fc — l)'th level {Rk-i{t)). How- 
ever, in the uncoupled region (with size Rjjit)), there 
is only the feeding from the boundary of Rcit). So in 
Ru{t), the cluster evolves with its own reaction dynam- 
ics. It implies that the number of particles in the fc'th 
level (iV/c(i)) increases in time with different exponents at 
the multicritical point, viz., Nc ~ t'' and Nu r-j fi k 
in the coupled and uncoupled area respectively. It is 
clear that this inhomogeneity of the total particle num- 
ber Nk{t) = Nc{t) + Nu{t) is a general property of uni- 
directionally coupled systems. 

The significance of this inhomogeneity is illustrated by 
noting that at the multicritical point, each level decays 
more slowly than its lower levels. So the coupling to lower 
levels is always broken after some time. To exploit the in- 
homogeneity of Ni-{t) for any fc > 1 we only consider the 
configurations in which the source (fc = 1) still survives. 
In these configurations, slave systems {k > 1) always 
survives with unit probability so the survival probabil- 
ity Ps{t) of the slave is just that of the source and we 
have 6k ^ cxk and 6k = 6\ for any fc > 1. We call these 
configurations as source ensemble and the average with 
this ensemble as source average. As shown below the 
source average exhibits more clearly the special proper- 
ties of unidirectionally coupled systems than the usual 
slave average and allows for the derivation of the scaling 
relations without any assumptions of scaling functions 




FIG. 1: The evolution of a cluster of the second level of 
CP-CP coupling. Black(gray outside black) corresponds to 
coupled(uncoupled) region. 

for observables. 

In a coupled region, Nc{t) is given by the relation, 
Nc{t) = pc{t) X Rc'^it) X P/=^(t), where pcit) and 
Ps'^^^it) are the particle density of fc'th level inside of 
the volume Rc'^{t) and the survival probability of the 
first level. Assuming the power-law decay of pc, Rc 
and P/=i as pc f-""", Rc = Rk-i t^^'"'-' and 
Ps''~^ ~ t~^^ , we obtain the first scaling relation 

rj'^k = d/zk-i - a'^k - Si . (5) 

In the uncoupled region, p^ ^ and Ru ~ t^/'^"''. 

Then similarly we find the second relation 

r?^, = d/z^'k - a^'k - Si . (6) 

From Rk = Ru + Rk-i ^ t^^^'', we have a condition, 
z'^k > Zk- The above two relations are generalized scal- 
ing relations for unidirectionally coupled systems taking 
the inhomogeneity into account. 

Our scaling relations are very different from Eqs. (3) 
and (4) because they are not scaling relations for quan- 
tities of the whole system. The number of total par- 
ticles Nk{t) and total density Pk{t) scale as Nk{t) = 
Nc{t) + Nuit) ~ t"" and pk{t) = pc{t) + pu{t) ~ i"""- 
Therefore rjk should be the larger one of rf^ f. and rf k^ 
while Qffe be the smaller one of a'-^ k and k- For 
z^' k = Zk-i, rf k is larger than rf ^ because oF k and 
oP k always satisfy oP k < ct'^fe- In that case, we have 
^fe = k and OLk = of k and we get r\k = d/ Zk+ ak+ Si 
with Q!fc ^ ^fc. This is of the same form as Eq. (4) for 

k = Zk-i but with ak p Sk- However the values of 
r]k is not the same in general because it is measured in 
different ways. 

Notice that the inhomogeneity of a cluster gives cor- 
rections, t~^'', t~^i and t^'^' with some negative co- 
efficients to the scaling with the asymptotic exponents 
ak, rjk and Zk- The quantities Ap, and A^ are de- 
fined as Ap = \a^k - oPk\, \ = \il'-''k - k\ and 
Aj, = d\l/z^k — l/-2fc-i| respectively. These correction 
exponents may be small so that they can cause a signif- 
icant long time drift of effective exponents measured in 



3 




FIG. 2: A schematic phase diagram of a two-level hierarchy. 
The region labeled by 'active(inactive)' corresponds to an ac- 
tive(inactive) phase of both levels. In the other regions, each 
level is in a different phase. 




FIG. 3: Plot of the effective exponent of the second level in 
CP-CP coupling at the multicritical point. The lines corre- 
spond to ric TjB and r]u (from top to bottom). 



simulations. This makes the precise estimation of expo- 
nents very difficult in unidirectionally coupled systems. 

We test the scaling relations Eq. (5) and (6) in a two- 
level hierarchy in one spatial dimension. Each system 
belongs to the DP or PC class. We call the first and the 
second level A and B respectively. Then for the second 
level we have the scaling relations 

rjc ^ 1/ ZA - ac - 5a , 
rju = l/zu ~ au - 5a , 

where we omit the index B from rj and a. We denote the 
coupling of two systems as Source-Slave coupling, where 
the source system feeds particles to slave one. 

There are four ways of coupling DP and PC system. 
For PC-PC coupling, the existence of a multicritical point 
depends on the way of coupling 0|. We will dis- 
cuss it in other place. For DP-DP coupling, we consider 
a two-level hierarchy of contact process (CP) 0|. In 
the first level, a particle A is spontaneously annihilated 
with probability p and it creates one A particle on one 
of the nearest neighboring (nn) sites with (1 — a){l —p). 
In the second level, the annihilation and branching oc- 
cur with p and (t(1 — p) respectively. In branching, if 
the target site is already occupied by an other parti- 
cle, the branching process is rejected. The probability 
a is introduced for the coupling. The criticality of each 
level is pc^ = Pc^ = 0.232674(4) at cr = for the 
first and cr = 1 for the second level \\^. Using the fact 
that the ratio of creation and annihilation probabilities, 
i? = (1 — pc)/pc should be same at criticality for any 
value of cr, we find that the critical line of each level 
in cr — p phase diagram is (Jc^ = 1 — pR/ {1 — p) and 
a^B ^ pji^i _ ^j^i^ ^ ^ 3.29785 (Fig. 2). The 
coupling dynamics A ^ A A- B with (t(1 — p) linearly 
and unidirectionally couple two CPs without feedback 
from the second to the first level. The multicritical point 
{o'MtPm) is the intersection point of the two critical lines, 
CTc^ and CTc^. We find ctm = 1/2 and pM = 0.13165. 
If the multicritical point is approached along the line 
A = {pm — p) — > with fixed a = aM, new critical 



behavior emerges p. 

By measuring rjc and rju of Eq. (7), we can pre- 
dict rjB and as defined as Nb = Nc + Njj t^^ and 
Pb — PC + Pu ^ t~°"^ ■ Using source average at A = 0, 
we measure effective exponents ric{t), rnjit), ac{t), r]B(t) 
and l/zu{t) up to 10^ Monte Carlo time steps. r]c{t) is 
defined as r]c{t) = \og{Nc{mt)/Nc{t))/logm and sim- 
ilarly for the others. Table shows the numerical results 
of source average, r^s is fairly underestimated due to the 
large correction with A^ = 0.08(1). From Eq. (7) we 
predict ac = 0.0725(100), au = 0.1525(125). Because 
of rjc > rjif and ac < au, we predict rjB = rjc = 0.40(1) 
and aB — ac = 0.0725(100). The prediction agree very 
with the numerical results. Fig. shows rjc (t), rju {t) a,\id 
rjB{t). Our results are also equivalent to the previous re- 
sult of the second level of DP-DP coupling measured in 
slave average, rjB = 0.39(2) and as = 0.075(10) :Sj. 

Next, we consider a two-level hierarchy of systems be- 
longing to DP and PC class. We only discuss PC-DP 
coupling here. 

We consider the CP and branching annihila ting ran- 
dom walks with two offsprings (BAW(2)) [ij [13 ■ 
BAW(2) belongs to PC class. In BAW(2), a particle, 
A hops to one of the nn sites with probability p and 
it creates two particles on two nn sites to the left or 
right direction with (1 — cr)(l — p). If two particles 
happen to be on a same site by branching or hopping, 
they annihilate each other instantaneously. The criti- 
cal point of BAW(2) at cr = is at p^"^ 0.5105(7) 
[l3| . The critical line is = 1 — pi?^/(l — p), where 
= (1 -p/)/p/ = 0.95886. In CP, annihilation and 
branching of a _B particle occur with p and cr(l — p). 
Criticality of CP at ct = 1 is same as before. The critical 
line is now (Tc^ = pR/{l —p)- The dynamics A ^ A-^ B 
with (t(1 —p) couple BAW and CP without feedback from 
CP to BAW. An A particle creates a B particle on the 
same site of the CP level only if the target site is empty. 
The multicritical point {(Jm,Pm) of BAW-CP coupling is 
= 0.77474 and pm = 0.19023. 

At the multicritical point, we estimate rjc = 0.14(2), 
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r]ij = 0.1875(25) and 2/zij = 1.26(1). Using random ini- 
tial conditions of pa{0) = Pb{0) = 1/2, we also measure 
the exponent /3 of p"^ in steady-states up to system size 
1.2 X lO'' and A = 6 X IQ-^. We estimate Pb = 0.28(1). 
The numerical results show that the CP still belongs to 
the DP class. 

With ZA = 1.750(5), (Sa = 0.285(2) for BAW(2) [13, 
Eq. (7) predicts ac = 0.146(20) and au = 0.1575(96). 
The value of ac and ajj agree with the DP value a = 
0.15964(6) ^ within numerical error. It means that ps 
uniformly decays but very slowly converges due to the 
coupling effect in the coupled region. Due to rjc < Vu, 
Eq. (7) predicts rjg — rju = 0.1875(25) which is not 
the value of DP class, 77 0.313686(8) i^]. We directly 
measure -qB = 0.18(1) which is somewhat underestimated 
but agrees well with the prediction. We conclude that 
even though CP is not changed by BAW(2), Eq. (7) is 
satisfied very well while Eq. (4) gives a wrong prediction 
for this coupling. On the other hand, with the DP values 
for -qs and z_b, Eq. (4) predicts as — 0.0338, which is 
completely different from the true DP value, a = 0.15964. 
Eq. (4) also gives a wrong prediction for BAW-CP-CP 
coupling in which BAW(2) is on the first level. The CP 
of the third level is affected by the second level. But 
the second level is not affected by the first. In slave 
average, exponents of the third level are same as those 
of the second level in CP-CP coupHng without BAW(2). 
These results cannot satisfy Eq. (4) because 5i is not a 
DP value. 

In order to interprete these observations we stress that 
the dynamic exponent 77 can depend on average meth- 
ods, but a does not because it is the ratio of [3 and 
which characterize ofF-critical behavior. Since away from 
the critical point the survival probability for source is 
unity and hence the source and the slave average coin- 
cide. Therefore the results of two average methods should 
give the same value of a. So it is natural that two values 
of a measured in two different ways are same in CP- 
CP and BAW-CP coupling. The method of average is a 
matter of choice but one obtains distinct scaling relations 
which are valid for any coupling. It means that Eq. (4) 
is not a relation for the slave average because it is not 
always satisfied by the results of slave average. Also Eq. 
(4) is not a general relation for unidirectionally coupled 
systems because it is a relation for quantities of the total 
system. As the relation for quantities of a total system 
is obtained by comparing Eq. (5) and (6), it can change 
according to the nature of the coupling. 

In summary, unidirectionally coupled systems show an 
inhomogeneity in the number of particles or density. A 
cluster in a higher hierarchy is divided into a coupled 
and an uncoupled region. In each region, dynamic quan- 
tities show different scaling behavior. In order to describe 
this property quantitatively we employ a special average 
method, the source average. Using scaling arguments we 
derive two generalized scaling relations which describes 



TABLE L Estimates of source average for the dynamic ex- 
ponents of the CP-CP coupling. 

nc Vu '2/zu ac Vb 

0.40(1) 0.32(1) 1.265(5) 0.075(5) 0.375(25) 



the inhomogeneity. The presence of 6 of the first level 
in the scaling relations arises naturally from the source 
average. Our scaling relations also explain the existence 
of large corrections to scaling of dynamic quantities of a 
total system which cause a long time drift of critical ex- 
ponents. Our scaling relations should be tested for other 
couplings between different universality classes and also 
in higher dimensions. The DP-PC coupling will be a 
possible candidate for future study. 
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ship Program of Korea Science & Engineering Founda- 
tion (KOSEF). 
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